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ABSTRACT 

A numerical double burst model of the chemical evolution of gas-rich dwarf galaxies 
has been developed. The model is fitted to a sample of N/O, O/H, Y and gas fraction 
observations, where N/O and O/H are the relative abundances by number of nitrogen 
to oxygen and oxygen to hydrogen, respectively Y is the abundance by mass of helium. 
Closed models as well as models including enriched outflow, ordinary outflow and 
ordinary outflow combined with inflow are considered. The bursts are assumed to be 
instantaneous but ordered in pairs to explain the scatter in N/O-O/H. The method 
of gas fraction fitting is revised, and it is found that it is very important to specify 
whether dwarf irregulars (dlrrs) or blue compact galaxies (BCGs) are considered. 
Effective enriched winds fail when fitting N/O, whereas closed models, models with 
ordinary winds or a combination of ordinary winds and inflow are all viable. 

Key words: methods:numerical - galaxies: irregular - galaxies:starburst - galax- 
ies: abundances - ISM:H n-regions 



1 INTRODUCTION 

The chemical evolution of dwarf irregular (dlrr) Q and blue 
compact emission-line galaxies (BCGs) is of particular inter- 
est because a substantial body of observational data is avail- 
able and some degree of simplicity exists because of the low 
level of 'metal' enrichment and absence of large abundance 
gradients. Furthermore, their wide range of intrinsic proper- 
ties makes them suitable objects for testing certain expecta- 
tions from stellar nucleosynthesis theory and the 'Simple' or 
other models of galactic chemical evolution, although at the 
same time there are complications associated with inflow 
of unprocessed material, outflow in homogeneous or selec- 
tive galactic winds and bursting (or 'gasping') modes of star 
formation. Chemical evolution models attempt to apply all 
these concepts to account for the distribution of different 
elements, notably helium, oxygen and nitrogen, in relation 
to star formation rates and gas fractions. Because many pa- 
rameters such as these last two are generally very poorly 
determined, the most convincing tests come from the com- 
parison of different elements with one another. 



Back in the 19 70s, S mith ( |1975| ), Peimbert fll97gp and 
Edmunds & Pagel (1978) noticed a contribution of primary 
nitrogen to the N/O ratio in Galactic and extragalactic H II 
regions with low oxygen abundance and Edmunds & Pagel 
attributed the existence of scatter in N/O at a given O/H to 



* In this paper a dwarf galaxy is defined to have absolute mag- 
nitude M B > -17. 



the existence of a time delay in primary nitrogen production 
by intermediate-mass stars, combined with differing effec- 
t ive a ges of the underlying stellar populations. Alloin et al. 
( 1979| ) also noted the primary nitrogen and attributed scat- 
ter in N/O to variations in the initial mass function (IMF), 
whereas Lequeux et al. (1979) in their classic study of he- 
lium, nitrogen and oxygen in irregular galaxies and BCGs 
confirmed the primary nitrogen likewise, but were not con- 
vinced that there was any real scatter in their data. The 
models of Alloin et al. and Lequeux et al. assumed evolu- 
tion to tak e plac e smoothly as a function of time; Matteucci 
& Chiosi (1983) were the first to incorporate into chemi- 
cal evolution models for these systems the idea of burst- 
ing modes o f star formation as prevously inferred by Searle 
& Sargent (1972) and Searle, Sargent & Bagnuolo (1973), 
and interpreted on the basis of th e SSPSF hypothesis by 
Gerola, Seiden & Schulman (1980). The basic pattern of a 
primary (constant N/O) pattern at low metallicities in H II 
regions changing over to a secondary pattern (N/O tx O/H) 
at higher ones has been confirmed in many more recent in- 
vestigations (e.g. Vila-Costas & Edmunds 1993; van Zee, 
Salzer & Haynes 1998). 

The chemical evolution of dlrrs and BCGs has been 
studi ed in many more recent investigations. Matteucci & 
Tosi (1985) found good fits to the data with a Salpeter IMF, 



inflow, homogeneous outflow, bursting star formation and a 
choice of third dredge-up parameters from Renzini & Voli 
(1981), attributing scatter in the N/O ratio to variations in 
M up , the upper limit to the masses of stars undergoing the 
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third dredge- up with hot-bottom burning. Garnett ( 199C] ) 
indicated schematically how the occurrence of bursts could 
in itself lead to variations in the N/O ratio just as a result of 
observi ng systems at different stages in the burst cycle. Pi- 
lyugin (1992; 1993) developed similar ideas in quantitative 
numerical models involving self- enrichment of H II regions 
and selective galactic winds as well as burst ing m odes of star 
formation, and Marconi, Matteucci & Tosi (1994) also devel- 



oped models with bursts and selective winds, while Carigi, 
Colin & Peimbert ( 1999 ) have investigated similar models, 
but prefer a 'bottom-heavy' IMF similar to one claimed in 
some globular clusters and giving rise to low true yields. 
However, part of the motivation for invoking selective winds 
was the apparen t exis tence of a large dY/dZ ratio suggested 



1 ot al a mi ), whioh no longc 



valid 



;otov 



1998), and the scatter in N/O also seems to have 
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been overestimated in those investigations. In this paper, 
therefore, we investigate the problem again, making use of 
more recent data and models of stellar nucleosynthesis and 
exploring in particular the role of mixing processes and of 
differing burst phases in leading to scatter in the N/O, O/H 
relation and the relationship between oxygen abundance and 
gas fraction. 

The structure of the article is as follows: Section ^ 
presents the adopted sample of abundance observations. Sec- 
tion |^ discusses the evolution of Hn-regions, wind-driven 
bubbles and supernova-driven supershells to investigate pos- 
sible mixing scenarios. This leads to the description of our 
double-bursting models in section ^J. The results of fitting 
the models to the sample of observations are presented in 
section [|, and discussed in section ^. 



2 THE OBSERVATIONAL DATA 

It has been stressed that the sample should consist mainly of 
BCGs, as these objects are extreme in both star formation 
rates and metallicity, thus being attractive for a bursting 
model. Another demand has been to include the most re- 
cent data only, in order to get the observations as accurate 
as possible. Further, a few damped Ly-alpha systems (DLA) 
have been included, mainly for the sake of comparison. The 
motivation for including them is found in the growing opin- 
ion that they represent dwarf galaxies in early (high z) stages 
of evolution. Thus, they may offer some insight into the ex- 
treme low-metallicity environment in the early evolution of 
present -day dwarf galaxies. 



Th 3 following sample is selected: BCGs and dlrrs from 
Pagel e t al. ( |1992| ), BCGs from Izotov, Thuan fc Lip ovetsk y 
( [19974 , the BCG SBS 0335-052 from Izo tov e t al. ( (t997b| ), 
the BCG IZwl8 from Izotov & Thuan ( |l998|). and finally 
four DLAs from Lu, Sargen t & Barlow (1998). From the 
sample of Pagel et al. ( |l992[ ) are excluded LMC, NGC5253 
and NGC5455/NGC5461 (the last two are Hn-regions in 
the nearby spiral M101), and objects included from one of 
the other sources. The total sample is shown in fig. jj] for 
N/O-O/H and fig. | for Y-O/H. The N/O abundances seem 
to have a constant level for 12+log(0/H) less than 8, im- 
plying the major part of the nitrogen to be produced as 
a primary element. For 12+log(0/H) higher than about 8, 
N/O is increasing as a function of metallicity, suggesting the 
major part of the nitrogen to be secondary. These trends are 
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Figure 1. The observational samp le used in this pape r. 
The symbols are: di amond=IZwl8 f|botov fc Thuan 199ij ), 
square=SBS0335-05 2 ( |lzotov et al. 1997b| ), filled circles are BCGs 



I Izotov et al. 1997a ) and open circles are BCGs a rid dlrrs fjpagcl| 



et al. 1992). Finally, the plusses are DLA-systems (Lu et al. 199£ ) 



Figure 2. The helium abundance sample. The symbols are as in 
fig. [J No DLAs are shown. 



wellkn own, see e.g. Garnett ( 199C ) , Vil a-Cost as & Edmunds 
( [l993| ), ! D ettini, L ipman & Hunstead ( |1995| ) and van Zee, 
Salzer & Haynes (|l998a| ). 



Further, the scatter in N/O seems to be significant and 
our models are based on this assumption. If we fit a straight 
line to our selected observations, ignoring the DLAs, we find 
the standard deviation in N/O to be cr=0.11, which is of the 
s ame order of magn itude as the observational uncertainty 
(Izotov et al. 1997e), but a real scatter is evident in other 
data from a wi der range of sources, see e.g. Kobulnicky & 
Skillman ( |l996| ). 

The DLA-systems pose severe problems when observ- 
ing N and O abundances. Usually the N-lines are occuring 
on top of underlying absorption, and O-lines are almost al- 
ways satura ted. Because of these difficulties, Lu et al. (Lu 



et al. 1998) found it useful to use Si or S instead of O, 



which makes sense as O/Si and O/S are found to be the 
same as solar in both the Galactic disc and halo as well as 
in Hn-regions in nearby galaxies. The question is whether 
the abundance ratios are equal to solar at the extreme low 
metallicity of DLAs. This is by now the largest uncertainty 
in this method. Thus, in fig. [l], the DLAs are represented by 
their N/Si and Si/H, assuming that these ratios are equiv- 
alent to N/O and O/H. Four systems have been selected 
from the Lu et al. sample namely those in front of QSO's 
0100+1300, 1331+1704, 1946+7658 and 2343+1232. For the 
rest of their sample, the quality of the observations restricts 
the abundance determination to be performed as higher or 
lower limits. Lu et al. did not make any corrections for dust. 
However, they noted that the depletion in dust is less than 
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0.4 dex. Though the observations of the four DLAs are in- 
cluded, they should not be taken too seriously and are not 
given much attention in this paper. 

The helium mass fraction shows a linear dependency 
on metallicity. This linearity has been used for extrapola- 
tion back to O/H=0 , giving the primordial He abundance. 
Izotov et al. (1997a) used their observations, included in 
the present sample as well (filled circles in fig. ^), to obtain 
0.243±0.003, using a linear fit with slope dY/dZ=1.7±0.9 
and assuming Z=20(O/H). Error bars are omitted in the 
figure to keep the clearness. However, the error bars are in- 
cluded in the figures presenting our results in section ^ and 
show no observational evidence for a scatter. Omitting SBS 
0749+568, represented by the isolated point above 0.26 in 
fig. ^| but including all other galaxies in the sample, a linear 
least-square fit gives dY/dZ=2.63±2.21, again using Z=20 
(O/H), and assuming a confidence interval of 95 per cent. 
This gives a primordial helium abunda nce y„ =0. 238=1= 0.004, 
which is consistent with Izotov et al. (1997a) within uncer- 
tainties. 



3 EJECTA DISPERSAL AND MIXING 

Many chemical evolution models assume a one-zone descrip- 
tion with instantaneous mixing. Assigning the term one-zone 
to a BCG seems to be quite a poor approximation, and as- 
suming the mixing to be instantaneous will always be doubt- 
full. However, if we choose to reject the assumptions, we are 
faced with the problem that no complete theory of mixing 
is available at present. The problem is that the gas dynam- 
ics following a starburst are very complex. The energy input 
into the ISM comes from photoionization, stellar winds from 
massive stars followed by SN type II. Different scenarios have 
been proposed for the dispersal and mixing processes. One of 
the propositions has been 'self-enrichment', referring to the 
suggestion that only the H n-region, surrounding the newly 
formed stellar cluster, is enriched with heavier elements, in 
particular oxygen. Thus, the observed abundances, using 
emission lines, may be considerably higher than they would 
be if observed in the neutral medi um. T he idea was originally 
proposed by Kunth & Sargent (1986), suggesting that the 
enrichment is confined to take place within the Stromgren- 
sphere. However, the hypothesis is controversial, and a lot 
of opposing arguments exist. Use of self- enrichment implies 
the assumption of almost instantaneous mixing within the 
Hn-regions. The problem is that the correlated energy in- 
put of SNe changes the physical conditions of Hn-regions 
dramatically. In fact, the density is decreased by a factor 
10 3 — 10 6 and the temperature is increased by a similar fac- 



tor behind the SN-driven shock front. The observed lines 
of single or double ionized N and O cannot arise from such 
extreme conditions. A future project could be X-ray abun- 
dance observations. If the observed emission lines cannot 
form within the superbubble/supershell (hereinafter a 'su- 
perbubble' is a wind-driven bubble, while a 'supershelP is a 
SN-driven shell), they may arise from a region outside the 
shock front. But according to detaile d numerical hydrody- 
namical models ( Tenorio-Tagle 1996 ), the ejecta will stay 
within the superbubble/supershell for a large part of their 
evolution, not enriching the surrounding medium. To test 
the viability of this statement, a few calculations are pre- 



sented below, comparing radii of superbubbles/supershells 
and radii of Hn-regions. Further, observations give evidence 
against the self-enrich ment hypothesis, see e.g. van Ze e et al . 
( |l998b| ) , Kobulnicky ( |l997| ), Kobu lnicky & Skillman ( |l997| ) 
and Kobulnicky & Skillman (|l998j). 



3.1 H II-region evolution 

In the following, the radii of Hn-regions are calculated 
for comparison with superbubble/supershell radii. The 
Stromgren-sphere is defined to be the sphere within which 
all ionizing photons are absorbed. Thus, setting the flux of 
ionizing photons equal to the number of recombinations, in- 
tegrated over the entire Stromgren-sphere gives: 



47T -i 

—R(tfa B n e n p € = N(t) 



(1) 



e is the filling factor, as is the reco mbination rate e qual 
to 2.59 x lO-^cmV 1 (T=10000K) jOsterbrock 1989| ). n v 
and n e are the number densities of protons and electrons, 
respectively, and N(t) is the flux of ionizing Ly c -photons. 
Assuming that all hydrogen inside of the Stromgren-sphere 



is fully ionized, n e 



The filling factor is an indicator of 



the uniformity of matter, having a value of 1 when the mat- 
ter is completely uniformly distribut ed. A typical va lue of 
the filling factor is of the order 0.01 (Kennicutt 1984). The 



electron density in H n-regions is typically n e ~ 100 cm 
(izotov et al. 1997a). When inserting these two values into 



eq. |lj, the corresponding rms density is n e — 10 cm , since 
e = 1 and n e = 10 cm -3 is equivalent t o the insertion of 
the observed values. According to Spitzer (1975, p. 251), the 
first phase in the existence of an H n-region is characterized 
by an almost static Stromgren-sphere. Henceforth, the ra- 
dius of this initial sphere will be referred to as 'the initial 
Stromgren radius', Rq. During the second phase, the sphere 
is expanding until the H n-region vanishes. If the flux of ion- 
i zing p hotons is assumed constant for the moment, Spitzer 
( 1973) found the relation 



R <'>=»'( i+ i^r 



(2) 



on the basis that the expansion velocity nearly equates the 
velocity of the associated shock, t is the time since the re- 
gion started to expand, virtually equal to the age of the 
burst, and Ch is the sound speed in the Hn-region , equal 
to 17kms~ 1 (T=10000K, 7 = 5/3) ( jOsterbrock 198S| ). Eq. § 
is not realistic, though, because the photon flux is decreasing 
in time. In the following this is taken into account, extend- 
ing the calculations by Spitzer. Differentiating eq. |l| with 
respect to time gives 



dN -iRT) 2 r? dR j_9rtB 3 „ d P u 



(3) 



where (3 = ^■rfa.Bf-, n e = r/pu, so r\ gives the relation be- 
tween the number density of electrons and the mass density. 
Below it is shown that 77 cancels in the calculations, so we 
will not worry about its value. Finally pu is the mass den- 
sity within the Hn-region. Consider a shell of ionized gas. 
Assuming uniform expansion, i.e. 



1 dr 
r~dt 



Vi 

R 



(4) 
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where r is the comoving radius of the shell and Vi is the 
velocity of the ionized gas just within the ionization front, 
assumed to have the same radius R as the Stromgren-sphere. 
The mass inside the comoving shell is conserved by defini- 
tion, so r 3 pn is constant leading to 



1 dr 
r~dl 



1 dp! 



3pn dt 

Inserting these two equations into eq. 



gives 



Vi 



1 dR 
2~dt 



'ILL 

dt 



&PR 2 P 2 n 



(5) 



(6) 



and 



— Vi, the velocity of the ionization front. The phys- 
ical conditions across the shock, formed ahead of the ioniza- 



tion front may be described by the jump condition 
pi V 2 = pi i + piiu 2 



(7) 



pi being the mass density of the surrounding neutral 
medium, V s is the velocity of the shock, pn is the pres- 
sure within the Hn-region and u s is the inward velocity of 
matter with respect to the shock. It has been assumed that 
the density and pressure between the shock and ionization 
fronts are the same as within the ionization front, and the 
pressure of the surrounding neutral medium has been ne- 
glected. Assume V s = Vi, so u a = m — Vi — Vi, where m is 
the inward velocity of matter with respect to the ionization 

front. Inserting this and the equation of state pn = 11 11 
with 7 = 5/3, pii = y/ (from eq. [j]) and pi = y/^rj 

(i.e. the density before expansion starts) into eq. [7] and ^j, 
one finally arrives at the second order equation 



NRl \ 
N R? J 



-1/2 



(dR\ 2 f (dR 
\dt J WR 2 P 2 u V dt 



/ dN_ 
I dt 



\6pR 2 Ph 



+ 



r 2 
^11 



= 



It is straightforward to solve this equation for 4§, obtain- 
ing a differential equation, which is solved numerically, us- 
ing timesteps equal to 1 Myr. For the calculations we used 
n e — 10 cm" 3 (/3 = l)as the initial density. When inserted 
into the second order equation instead of pn, the factor r\ 
cancels remembering that it is included in f3 as well. For 
every timestep, the density ne.{t) is calculated from eq. [ij 
The fluxes of ionizing p hoton s are taken from the models by 
Stasinska & Leitherer (1996), giving the fluxes from 1 Myr 
to 10 Myr after the burst in intervals of 1 Myr. An instan- 
taneous starburst of mass 1O 6 M was assumed, comparable 
to those used in our model. Further, a Salpeter IMF with 
the same upper mass limit is used, only demanding a sim- 
ple scaling from their lower mass IMq to our lower mass, 
either 0.1 or O.OIMq. If mr, = 0.1 the normalization con- 
Integration of the IMF yields 



stant is 0.17, see section 4.1 



a mass fraction of 0.39 above 1 M@, so the adopted fluxes 
are multiplied by 0.39. Correspondingly, using tul = 0.01 
the normalization constant is 0.07, giving a mass fraction of 
0.17 above 1 Mq. The decrease of continuum photon flux 
calculated by Stasinska & Leitherer using their parameter 
values is shown in figure ^[ 

Ro is found from eq. ^ by inserting the initial flux of 
photons and the initial density n e = 10 cm -3 . The initial 




Time since burst (iv'yr) 



Figure 3. The decrease in ionizing photon flux counted from a 
starburst event. 



flux, No = N(t — 0) is the flux just before the Hn-region 
starts to expand. No information exist on the photon flux 
before 1 Myr after the burst. Thus, No = A^(lMyr) is as- 
sumed so Ro is equal to R(l Myr) found from eq. |l[ The 
calculations are performed for both of the IMF's and the 
results are shown in table 



3.2 Supershell evolution 

In the case of a single supernova event, it is rather simple to 
follow the four main evolutionary phases of its remnant ana- 
lytically. First, the ejecta and the leading shock move almost 
undecelerated. Second, the velocity of the shock decreases, 
as a shell of swept-up ISM builds up and grows larger. In 
this phase, radiation from the hot cavity is not important, 
and the total energy within the shock is conserved. Third, 
radiative cooling of all matter within the shock become im- 
portant. The shell is now driven by the pressure still existing 
in the cavity. Finally the PdV work from the cavity vanishes, 
so that the shell is now only being driven by the momentum. 
At some moment, the velocity of the shell has dropped to a 
value comparable to the rms velocity of the ISM, thus being 
non-detectable. 

When one deals with the correlated effect of many su- 
pernovae, which is the case in starburst galaxies where the 
number of massive stars is large, additional considerations 
have to be made, since supernovae keep feeding energy to 
the cavity for about 50 Myr. The idea with the following 
calculations is to get a feeling of the timescales and radii 
of supers hells. The analytical expressions from McCray & 
Kafatos (1987) are adopted to do the simple calculations, 
and the equations below are taken from their article, un- 
less another reference is given. Strictly speaking, the ana- 
lytical expressions may be used only if the ambient density 
is constant. However, the supershell first expands into the 
wind-driven bubble environment, then into a region with 
density depending on the Hn-region expansion, and finally 
into a non-expanding region outside the H 11-region. Hence, 
the evolution of the supershell is only followed while it is 
smaller than the H n-region, and the density within the su- 
perbubble is assumed to be the same as outside. This last 
assumption is not realistic, as the passage of the wind-blown 
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shock sweeps a large fraction of the gas into a superbubble. 
Hence, we may expect a more accurate model to show that 
the supershell reaches the superbubble at an earlier time. 
The gas gathered by the wind-driven shock cools to form 
the superbubble after relatively short time 



Table 1 . H Il-regions and supershells 



i fr, o 1 n 4 \ -0.71 r 0.29 

tb = (2.3 x 10 yr)n e L 38 



(9) 

( Mac Low fc McCray 198^ ), and the radius of the superbub- 
ble is 



R s = 269 pc (L 3 8/n e ) 1/5 4 /5 



(10) 



where L38 is the characteristic wind luminosity in units of 
10 38 ergss _1 and is the elapsed time since the starburst in 
units of 10 7 yr. Strictly speaking, the time is counted from 
the start of superbubble expansion, but since O-stars start 
their wind-phase soon after the burst, t~j is counted from 
the burst. n e is the matter density within the Hn-region, 
assumed to be homogeneo us. T he wind luminosity is found 
in Leitherer & Heckman (1995) from their fig. 55 (instan- 



taneous starburst, metallicity is one tenth of solar, mass of 
burst is 10 6 M o , Salpeter IMF with upper limit 100 M© and 
lower limit IMq). At 10 6 yr the wind luminosity is read- 
off to be 10 39 ergss _1 . Doing the same IMF scaling as in 
the photon flux discussion above, L38 = 3.9 is obtained if 
m L = O.IA/q and L 3 8 = 1.7 if m L = 0.01M©. The wind 
phase is important for as long as O-stars exist (« 5 Myr). 
However, after about 3 Myr the first supernovae explode 
and soon dominate the total energy output. The superbub- 
ble radius has been calculated 1, 2 and 3 Myr after the burst 
using eq. [l(]. The density typical for an Hn-region during 
its early evolution is 10 cm -3 , as discussed above. 

After 3 Myr, the first SN appear. Hence, the expansion 
of the supershell is set to start at t=3 Myr. Strictly speak- 
ing, the energy output from supernovae appears as discrete 
events. However, it can be treated as continuous, as long as 
the interval between explosions is sufficiently short, at most 
in the order of 10 5 yr (Tomisaka, Habe & Ikeuchi 1981). If 
the total mass of the burst is 1O 6 M0, the number of stars 
with mass greater than or equal to 8Mq (lower lim it of a 
SNII progenitor, according to Woosley & Weaver (1986)) 



may be calculated using our IMF. The result using the 0.1 
Mq lower mass limit is 8145 stars, whereas it is 3452 stars 
using ttil — O.OIMq. The event lasts for about 5 x 10 7 yr 
so the mean interval between explosions is in the order of 
6 x 10 3 yr and 1 x 10 4 yr, respectively. Thus, it is meaning- 
ful to assume continuous energy injection. Using this, eq. hd 
may be used with just one replacement, namely the insertion 
of the mean supernova power instead of the wind luminosity. 
The mean power is calculated by dividing the total SN en- 
ergy iVi-Esi x 10 51 ergs with the total duration of the event. 
Nt is the number of stars with mass greater than 8 Mq and 
E51 is the energy in units of 10 ergs. Here, 0.4 is ad opted 
as a representative value, see Woosley & Weaver (1986, their 
table 1). The energy ejection is assumed not just to be con- 
tinuous, but also constant in time. This is probably quite a 
good approximation, see Leitherer & Heckman (1995). The 
considerations lead to 



R 3 = 97 pc (N t E 51 /n e ) 1/5 t 3 7 /5 

V s = 5.7 km s _1 (./V*£ 5 i/n e ) 1/5 t~ 2/5 



(11) 



t(Myr) 
1 
2 
3 
1 
5 
6 
7 



%n (P c ) 
209 (158) 
222 (172) 
233 (184) 
239 (189) 
245 (196) 
250 (201) 
255 (206) 



R a (pc) 
56 (47) 
85 (72) 
108 (92) 

123 (103) 
186 (157) 
237 (200) 
282 (238) 



n e (cm 3 ) 
10.0 (10.0) 
8.8 (8.5) 
7.5 (7.1) 

4.4 (4.1) 

3.5 (3.3) 
2.5 (2.3) 
1.8 (1.7) 



The radii and densities of H Il-rcgions compared to the radii of 
superbubbles/supershells at different times after the burst cal- 
culated with the rrt£=0.1AfQ cutoff. Values of shell radii typed 
in bold, are the radii of the superbubbles. For t>3 Myr, the 
calculated shell radii are for the supershells. The values for 
m£=O.OlA'/0 are given in brackets. The size of the superbub- 
ble at 4 Myr after the burst is roughly the same as the size of 
the supershell, so the supershell reaches the superbubble about 1 
Myr after the appearance of the first SNe. 

The time of shell formation is given by an expression similar 
to eq. ^, obtained by replacing the wind luminosity with the 
mean supernova power. 

The density of the Hn-region is assumed constant in 
time and equal to 3 cm -3 , a reasonable value during late 
stage Hn-region evolution (see table |l|). Eqs. [ll] are used in 
the adiabatic phase only (no cooling), since they implicitly 
assume that the loss of energy within the shell is negligible. 
Cooling, and hence energy loss, from the interior becomes 
important at a time 

i c = 4x 10 6 yr C^'iN+BkifV- 7 (12) 



At this time, the shell has a radius 
R c = 50 pc e-°- 9 (JV*£?Bi) -V°- 8 



(13) 



found by inserting eq. [li] into eq. [n]. £ is the metallicity in 
units of the solar metallicity. For the objects in question, 
£ = 0.1 is a representative value. Even if a small starburst 
mass of 10 Mq is assumed and the density of the surrounding 
medium is 3 cm -3 , it turns out that the supershell reaches 
the end of the adiabatic phase at a time 2.1 x 10 7 yr if 
rr>j,=0.1 is used and 1.6 x 10 7 yr if m,L=0.01. Hence, the 
interior of the supershell starts to cool after the H 11-region 
has vanished. This tells us that we may use eqs. to cal- 
culate the size of the supershell, as we are only interested 
in the evolution during H 11-region existence. The results are 
given in table |l| 

The measured emission lines may indeed originate in 
the part of the H 11- region that is still outside the superbub- 
ble, since the size of the H 11-region always exceeds the size 
of the superbubble as seen in table jjj. It also exceeds the 
size of the supershell until 6 Myr after the burst, though 
only significantly within the first 5 Myr. The above calcula- 
tions do not prove that the emission lines arise outside the 
superbubble/supershell - they just confirm that it is a pos- 
sibility. If emission lines really are measured this way, they 
are not affected by enrichment from the present burst, thus 
opposing self-enrichment. 

3.3 Propagating star formation 

Propagating star formation means that the energy deposited 
in the ISM by an evolving star formation event, is initiating 
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new star formation. As a supershell grows, it will become 
Rayleigh- Taylor unstable a nd fra gments. The situation has 
been treated by Elmegreen ( 1994 ) , who found an expression 
for the timescale of cloud-collapse in the supershell: 



= 103 



V cm _J / 



-1/2 



Myr 



(14) 



where M is the expansion speed of the shell divided by the 
rms velocity disper sion in the shell . For an adiabatic shell, 
M is equal to 1.8 (Elmegreen 1994). If n = 10 cm -3 , this 



gives the result that star formation is expected to start not 
earlier than ~ 24 Myr after shell formation. On the other 
hand, if no = 3 cm" 3 as used in the calculations above for a 
typical density in the late stages of H n-region evolution, the 
timescale is more like 44 Myr. The mass of the supershell is 
given as 



M = po V S! 



1.3nom p -^-R 3 



(15) 



assuming all of the ISM originally within R to be incor- 
porated in the shell, po is the mass density of the ambient 
medium, V ss is the volume occupied by the hot phase within 
the supershell, m p is the proton mass and 1.3 is the approx- 
imate mass per particle in units of the proton mass. The 
radius is of the order a few 100 pc as shown above, so eq. 
|l5| gives a shell mass of the order 10 6 Mq. Hence the in- 
duced star formation is of the same order of magnitude as 
the original central burst. 

Observational indications on the existence of star form- 
ing supershells are numerous. One recent exam ple is NGC 
2537 included in the sample of Martin ( 1998] ). In H a , it 
shows a very clear spherical distribution of starformation 
s ites. As a curiosity, it may be mentioned that Mori et al. 
(1997) calculated the effects of star formation in an expand- 



ing supershell, using a 3-D hydrodynamical code including 
a dark matter halo. The results were in remarkable accor- 
dance with available observations of dE's, such as exponen- 
tial surface brightness profile, positive metallicity gradient 
and inverse color gradients. 

If the ejected metals from the central burst mix with 
the material of the supershell, the Hn-regions of the in- 
duced burst will show abundances differing from the central 
H n-region. Although there have been some doubt whether 



the me tals are allowed to mix into the supershell ( Tenorio- 



Tagle 1396), and if so when this will happen, we have been 
inspired by this possibility to let our model have a double 
burst nature. 



4 THE MODEL 

We have designed a model to fulfil certain demands, namely 
that the observed trend of constant N/O for low metallicities 
and increasing N/O for higher metallicities should be repro- 
duced and also that the observed scatter in N/O should be 
explained. Fixing all the parameters by fitting the N/O-O/H 
evolution, the model should also be able to explain the ob- 
served helium mass fraction as a function of O/H and O/H 
as a function of gas fraction. 

The included elements are H,He,C,N,0. The produc- 
tion of N is still a hot topic, since the degree of primary 
production at various metallicities is unclear. Further, N 



is produced mainly by intermediate mass stars while O is 
produced by massive stars, making it necessary to include 
stellar lifetimes. However, the inclusion of both N and O 
provides the opportunity of constraining stellar parameters 
and mixing scenarios. 

To account for the scatter in N/O, it is suggested that 
the bursts are instantaneous and ordered in pairs, hence us- 
ing the delay between the ejection of N an d O. The principle 
of time delay was used by Garnett ( 1990] ), though employ- 
ing single bursts only. It is doubtful whether single bursts 
produce scatter, since at least some abundance observations 
have to be done at the time of N-release. However, this is 
between two bursts, where no giant Hn-regions are present, 
providing no possibility for abundance observations by emis- 
sion lines. Single bursts may produce scatter only if IMF 
parameters or enriched wind efficiencies vary from galaxy to 
galaxy. However, there is no indication of IMF variations, 
and the existence of enriched winds, on the whole, is ques- 
tionable, and may be important in extreme low-mass galax- 
ies only - see below. Double bursts produce the scatter quite 
naturally: It was found reasonable above that a localized 
burst results in another burst, this time in the expanding 
supershell, surrounding the original burst. The timescale for 
star formation in the supershell is found to be J; 2.4 x 10 7 
yr, comparable to the timescale for O-ejection, but shorter 
than the timescale for N-ejection, hence producing the de- 
sired scatter. 

In our model we assume that the very first burst is a 
single one, all other bursts appear in pairs. The interburst 
period between the two bursts of a pair is tuned to give max- 
imum scatter, found below to be 30 Myr. The time between 
two pairs is set to 1 Gyr. 



4.1 The IMF 

The IMF used throughout this paper is the single-power 
Salpeter-IMF: 



4>(m) = </> m (16) 

f alpeter 1955| ). <j>o is the normalization constant, found by 
" u m<p{rn)dm = 1. Using a lower mass cutoff tul = 0.1 Mq 
and an upper mass limit mjy = lOOM©, (f>o is equal to 0.17. 
These limits are standard values, very often used in the lit- 
erature. A Salpeter IMF still fits the observations (above 
about lA/(7)) quite we ll despite of its age and simple appear- 
ance (Leitherer 1998). 



The lower-mass cutoff is an important parameter be- 
cause the yield from a generation of stars is a function of 
the adopted lower mass. The reason is that stars less mas- 
sive than ~1 Mq do not eject metals, thus locking-up all the 
material from which they are formed. Hence, lowering the 
lower-mass cutoff implies a lower recycling fraction, giving a 
lower yield. As will be apparent from our results, it becomes 
attractive to invoke another value of tul namely 0.01 Mq. 
In this case the normalization constant is found to be 0.07. 
Thus, the yield will be a factor 2-3 lower when using the 
lower cutoff. 

Throughout this paper, a metallicity invariant Salpeter 
IMF has been used, since there are at present se veral obser- 
vational indications of an abundance invariance ( Wyse 199£ ; 
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Massey 1998). The number of stars with mass [mi ,7712] is 



given by 



N([m 1 ,m 2 ]) = 4>oM burst I m 2,35 dm (17) 

J mi 

Am 

mi = mj — 



m-2 



2 
Am 



where M bur3 t is the mass of a burst, i.e. the mass of gas 
turned into stars. The stellar mass grid may be as fine as one 
wishes. For use in the present model, a grid of 60 different 
stellar masses has been adopted, since this was found to give 
sufficient mass resolution for the present purpose. 

The mass of each burst may be constrained by the 
available observations of SFRs. The SFR for BCGs is typ- 
ically 0.1 - 1 A/ Q yr~\ For the mean SFR of a double 
burst to be within this range, the masses of the bursts 
may be estimated from M burst (l) + M burs t(2)/ti b ~ 0.1 — 1, 
where ta is the time between the two bursts of a pair. If 
t ib = 3 x 10 7 yr, and assuming M burst {l) ~ M burat (2), this 
gives M burs t ~ 1 .5 — 15 x 10 6 Mq, in agreement with Mar- 



lowe et al. (1995, their table 7) 



4.2 The equations 

The calculation of abundances is carried out just before ev- 
ery burst. The reason for doing this is simply that we wish 
to compare with the observed abundances, being measurable 
in H n-regions only. The adopted equations are described in 
t he fo llowing. The notation is similar to the one in Pilyugin 
(1993). The first step is to calculate the masses of each el- 
ement present in the gas phase just before every burst. For 
the j th burst: 



(18) 



Mi(tj) = Mi(tj-i) - AMj-iXi(tj-i) 

3-1 

k=i 
j'-i 

+ Y AM fc (Qi(T 3 -,fc) - Qifa-i,*)) 

tj refers to the time just before the j th burst. Xi(t k ) is the 
abundance of element i just before the fc* burst. Tj t h is de- 
fined as tj — tj,, so it is the time elapsed since the k th burst. 
Thus the element yields from stars with lifetimes shorter 
than Tj t k has to be included in the two summations. Q(rj,k) 



the (j — l) th and j th burst. The second summation term is 
corresponding to the first one, except for the fact that the 
mass of element i was newly synthesized. 

AWi(tj-i) is the mass of element i, leaving the galaxy 
as a result of starburst j — 1. This term consists of two 
parts: The part belonging to the ordinary wind, and the one 
corresponding to the enriched wind. Thus, the wind term 
may be written as 

AWi(tj) = WismAMjX^j) (19) 
+ W S N(Qi(t w ) + Q{t w )X t {tj))AMj 



(Pilyugin 1993). Wism and Wsn are the efficiencies of the 
ordinary and enriched winds, respectively. Both efficiencies 
are zero if the model is closed. The physics behind enriched 
winds is based on the principle that the supershell following 
a burst, breaks up, allowing the hot ejecta to blow out and 
escape from the galaxy. It is assumed that the wind is caused 
by supernovae type II, so the end of the wind phase t w is the 
lifetime of the least massive star exploding as a SNII, which 
is set to 8M0, consistent with the mass adopted for our 
supershell calculations. Hence, for stars with masses above 
this limit, a fraction Wsn of the ejected mass of element i 
is leaving the galaxy. For stars less massive, Wsn = 0. It 
is important to notice, in accordance with the physics in- 
volved that this efficiency factor is the same for all elements 
considered, but since oxygen is dominating the ejecta from 
SNII, the winds will be enhanced in oxygen. 

The other possibility is ordinary galactic winds, arising 
as a consequence of a general heating of the ISM, causing 
a fraction to leave the galaxy, i.e. the composition of the 
gas leaving the galaxy is the same as the composition of 
the ISM. The first term of eq. [19j is the instantaneous-burst 



representation of Hartwick-outflow (Hartwick 1976), giving 
a mass-loss proportional to the burst mass. By multiplying 
the mass of the wind with the fraction of the i th element in 
the ISM, taken just before the burst that is responsible for 
the wind, one obtains the mass removed of element i due 
to an ordinary wind. Wism,Wsn and the burst masses are 
treated as free parameters, though one restriction is made, 
namely that the yields of the first burst should be sufficient 
to place the second burst approximately at the abundances 
of IZwl8. 

AW 1 "' { T j,j-i)X' n f is the increase in mass since the 
last burst of element i in the ISM due to inflow of gas. The 
composition of the infalling gas, given by X i , is taken to 
be primordial. The inflow rate is given by 



M(t) 



M 



10 c -*/ T 

T inf 



(20) 



is the mass fraction of gas ejected from the k th generation ([Lacey fc Fall 1985 ; Sommer-Larsen & Antonuccio-Delogu 



of stars, just before starburst j. The composition is left 
unchanged by stellar nucleosynthesis, thus being the same 
as in the gas from which the stars were formed. Qi{rj t k) 
is correspondingly the mass fraction ejected of element i, 
but newly synthesized. It is clear from this notation that 
Q{ T j,j) — Qi( T j,i) ~ 0, since stellar lifetimes are finite. The 
first term on the right side is the mass of element i just 



1995 J. Mq is the total mass accreted for t 3> r in f and 
Ti n: f is the accretion timescale. The mass accreted at time 
t is found by integrating eq. [m] from to t finding 
M(t) = M (l - e~ t/T ">f), finally giving 



AW Inf {r 3 . 



before the (j 



burst. The second term is the mass of 



= M (e-^- l/T '"i - 



-uh 



(21) 
(22) 



element i that has been turned into stars at burst j — 1. The 
first summation term is the mass of element i ejected, with- 
out being changed, by the k th burst in the period between 



This is unfortunately giving two free parameters further, 
namely r in f and Mq. In models not including inflow 
AW Ini (T j:j -i)=0. The mass of gas just before starburst 
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number j is 

M B {tj) = Mgfa-ij-AMj-i-AWitj-i) (23) 

+ AW /n '(T«-i) 
j-i 

+ ^ AMk(Q(Tj t k) — Q(Tj-l,k)) 

fe=l 

Finally the mass of the entire galaxy is 

M dw (tj) = M^fe-r) + Af'^fe-i) - AWfe-i) (24) 

For the closed models Mdwitj) = Md tt (0). 

Qi and Q are adopt ed f rom models of stellar evolution 
as described in section 4.3. However, Qn(t) is calculated 
from 



E 



M*) = o 



(25) 



= H,He,C,N,0 



assuming all newly consumed H to go into the production 
of He, C, N and O. This simplification is realistic because 
the involved elements are by far the most important. 

From the above equations, it is now possible to calculate 
the abundances by mass in the ISM, just before a new burst: 

Mi{t 3 ) 



(26) 



However, the observed abundances are not given by mass, 
but by number. The abundance by number of, say 0(tj) 
relative to H(tj) is 



M Q {tj 



(27) 



M H {tj) A 

where A represent the atomic masses and M is the gas phase 
masses calculated by eq. |l8[ 

In the calculations it is assumed that the ejecta mix into 
the entire interstellar medium before the new burst appear. 
It should be noted that this assumption is doubtful for the 
short interburst period between the two bursts of a pair. The 
fact that the metallicity of second generation stars is slightly 
higher, is included in the sense that the adopted yields are 
metallicity dependent. The calculations are terminated if the 
amount of gas is insufficient for a new burst, or the age of 
the dwarf galaxy is more than 15 Gyr. 

The initial conditions to be put into the equations are 
the following. The very first burst (the single burst) takes 
place at primordial gas composition, i.e. Z=0, Xf/ e =0.243 
jlzotov et al. 1997a| ), X//=0.757, X c =X N =X o =0. 

Before the first burst, the galaxy consists of gas only. 
Hence, the total mass of the galaxy is equal to the mass 
of gas, in all calculations adopted to be 10 8 Mq. However, 
the situation is a little different for the inflow models. In this 
case, the initial mass of a galaxy is adopted to be zero. Grad- 
ually, the galaxy increases its mass by inflow of primordial 
gas, until the gas mass reaches a limit, sufficient to start 
the first burst. To be consistent, this mass limit has been 
adopted to be 10 8 Mq. So the only difference from models 
not including inflow is that the first burst is delayed by the 
time it takes to build up a sufficiently large gas cloud. 



4.3 The adopted yields and their implementation 

Since the present models follow the evolution from the very 
first burst at Z=0 to about solar metallicity, metallicity- 



dependent yields have been adopted. The Renzini & Voli 
(1981) yields for low- and intermediate mass stars are the 



most widely used for the purpose of chemical evolution mod- 
elling. Their strength is that they are giving the yields for 
several choices of convection and wind parameters. Further- 
more they give the primary yield of nitrogen (and 13 C) sepa- 
rately, which is very important for our purpose. Indeed, one 
of the major problems to solve for low metallicity objects 
is the extent of primary produced N, compared to the sec- 
ondary production. Marigo, Bressan & Chiosi (1998) made 
also a separation between primary and secondary produc- 
tion in the mass interval 4-5M© . The question whether stars 
more massive than 5Mq produce primary N is controversial. 
The question is whether it is necessary at all to introduce 
such a source for primary N production. This is one of the 
questions we address i n this paper. The yields by van den 
Hoek & Groenewegen (1997), also for low- and intermedi- 
ate mass stars, do not include separate results for primary 
nitrogen, so the yields given are a mixture of primary and 
secondary components. This is the reason, why we have not 
adopted their yields for more than just a reference, see sec- 
tion n 

For use in the models, yields ar e selected in a consistent 
way, i.e. yields from Geneva tracks ( Maeder 1992 ) and yields 



from Padova tracks (Marigo, Bressan & Chiosi 1996; Marigo 
et al. 1998; Portinari, Chiosi & Bressan 1997) are kept in 
separate sets of yields. The adopted yields are organized as 
shown in table Note that not all yields in sets 1, 2a and 
2b are deduced directly from the Geneva tracks, as va n den 
Hoek & Groenewegen fll997|), Renzini & Voli (Il98l|) and 



Woosley & Weaver (1995) use different indirect methods. 



The problem is to calculate the Q and Qi terms just 
before every burst, given the IMF and the ages of the pre- 
vious bursts from which the stars eject the elements. The 
mass ejected of an element i, at a time r after a burst is 
calculated as a sum of contributions from each star down to 
a stellar mass m T , corresponding to a lifetime r. The life- 
time of each star is found by linear Lagrange interpolations 
in both mass and metallicity between the values given by 
the stellar tracks. The following expression is used for cal- 
culating Qi 



Qi 



E 



Himj) N(mj) 



(28) 



where Nfjrii) is the number of stars with mass rrii, found 
with the adopted IMF, and qi(rrii) is the mass ejected of el- 
ement i from stars having this mass. The qiS axe found from 
the references in table H in two steps. First, a linear inter- 
polation is performed between the values for the available 
metallicities. If the metallicity is lower than the metallici- 
ties, for which the yields are known, primary and secondary 
components of N are treated differently. The yields for pri- 
mary elements are taken to be equal to the yields at the 
lowest metallicity with known yields, except for set no. 3, in 
which the primary N is obtained by extrapolating linearly 
from the two metallicities with known yields. For this set, an 
extrapolation should be more correct than just assuming the 
yield to be constant below the range of metallicities covered 
by the yield sources, since the range does not extend below 
Z=0.008 for the critical mass interval 4-5 Mq. The price 
to pay is an extrapolation reaching far beyond the covered 
metallicity range. Note that stars with masses less than 4 
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Table 2. The organization of the adopted yields. 



set no. 
1 



Mass range (Mq) 
9-100 
9-40 
0.8-8 



Metallicitics 
0.001,0.02 
0.002,0.02 

0.001,0.004, 
0.008,0.02 



Elements 
He,C,0,m rem 
Nsec 
He,C,N, 

,77~l rem 



Reference 
M92 (large mass-loss) 
WW95 

HG97 



2a 



9-100 
9-40 
1-8 



0.001,0.02 
0.002,0.02 
0.004,0.02 



He,C,0,m r em 

Nsec 
He,C,Nprim, 
Nsec,0,m rem 



M92 (large mass-loss) 
WW95 

RV81 



2b 



As 2a, but small mass-loss for stars with M > 9M ( 



© 



6-100 
4-5 

1.046-3 



0.0004,0.004, 
0.008,0.02 
0.008,0.02 



0.008,0.02 



He,C,Nsec,0,m rem 

He,Csec,Cprim, 
Nsec,Nprim,Osec, 

Oprim 
He,C,Nsec,0,m rem 



P97 



M98 
M96 



Sources: M92: flMaeder 1992| , His t ables 4.5 and 61. WW95:(|\\ooslcv fc Weaver 1994 thcir 
tables 5A,5B,10A and 10B), HG97:( |van den Hock fc Grocncwcgeii 1997| , thcir tables 13,17,35 
and 38. Tables 13 and 17 are using a wind parameter r?=4 for Z=0.008 an d Z=0.02. For 
Z=0 001 (table 38) and Z=0.004 (table 35), r) is 1 and 2 respectively), RV81:(Rcnzinifc Vol: 



1981 



their tables 3a,3d,3h and 3i. For all tables, ??=0.333. For tables 3a (Z=0.02, Y=0.28) 



and 3h (Z =0.004. Y=0,232), a i s equal to 0. In the two other tables, a is 1.5. an d mjjBB 



Mq), P97:( |Portinari et al. 1994 th eir tables 7 and 10) M98:( |Marigo et al. 1998 
4 and 5, m HB B=5 Mq) and M96:( [viarigo et al. 1996 , their tables 4 and 5) 



their tables 



Mq or larger than the upper mass for hot bottom burning 
(free parameter in the model for set 2a and 2b, and 5 Mq 
for set 3) are only producing secondary N. The secondary 
N yield is interpolated between the lowest metallicity with 
known yields and Z—0, since the yield of secondary N is at 
Z=0, according to the definition of secondary production. 
For H, He, C and O the yields interpolated/extrapolated 
are the sum of primary + secondary yields. This is without 
influence on the results. 

Second, the yields are interpolated linearly with respect 
to initial stellar mass, finally giving the qiS. For 2a and 2b, 
the mixing length parameter a, defined as the convection 
mixing length divided by the pressure scale height, is a free 
parameter. We know the yi elds for two different values of a 
from Renzini & Voli (1981), namely a=0 and a=1.5. Thus 
when implementing the yields, the qiS are found in three 
steps. Before the two steps described above are carried out, 
the yields given in the mass interval 4 Mq to the upper mass 
of HBB, rriHBB, variable between 5 and 8 Mq, are linearly 
interpolated with respect to a. 

The implementation of the mass ejected unprocessed is 
very similar. The expression for the total mass ejected since 
the starburst is written as 



Table 3. The true yields and return fractions for the sets in use. 



E 



(nii - m rem )N(mi) 



(29) 



m rem is the mass of the stellar remnant (white dwarf, neu- 
tron star or black hole) . Again the value of m rem is extracted 
from the sources by linear interpolation, first with respect 
to metallicity, then between initial stellar masses. 





Z 


PHe 


PO 


R 


Set 2a 


0.001 


0.032 (0.011) 


0.020 (0.007) 


0.274 (0.116) 




0.02 


0.031 (0.011) 


0.008 (0.003) 


0.283 (0.120) 


Set 2b 


0.001 


0.032 (0.011) 


0.020 (0.007) 


0.274 (0.116) 




0.02 


0.024 (0.008) 


0.014 (0.005) 


0.282 (0.119) 


Set 3 


0.001 


0.076 (0.026) 


0.014 (0.005) 


0.279 (0.118) 




0.02 


0.096 (0.032) 


0.014 (0.005) 


0.299 (0.127) 



The values given are for the high low-mass cutoff. Numbers in 
brackets are for the low cutoff. Note that the He yields for set 
3 are higher than for the other sets by a factor of about 3. The 
O yields are generally quite robust, though the yields of set 2a 
at solar metallicity (Z=0.02) are lower by a factor of about 2, 
compared to the other O yields at this metallicity. R is the return 
fraction. 

It is very useful to calculate the true yields of selected 
elements because they allow for an easy comparison be- 
tween the nucleosynthetic outcome of different stellar mod- 
els. The true yield of an element is defined to be the mass 
of an element ejected from a star generation, divided by the 
fraction locked up in stellar remnants and low-mass stars. 
The total mass ejected of an element, e.g. O is given by 
qo (rrii) N (mi) , where qo(jm) is found from the sources 
following the proces described above. The lock-up fraction 
is equal to 1— R, where R is the return fraction, i.e. the mass 
ejected in units of the total mass of the burst Q/Mf, urst . The 
resulting true yields are given in table pi for He and O, as 
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Table 4. The stellar lifetimes for the two stellar tracks. 



M 


Geneva 




Padova 


O 


M 


Zj — U.UU1 


17 ri no 

Zi — U.Uz 


Z=0.0004 


ry r\ rio 

Zj — U.UZ 


1 


100 


3.21E+06 


3.09E+06 


3.38E- 


-06 


3.39E+06 


1 


60 


4.08E+06 


3.89E+06 


4.19E^ 


-06 


4.12E+06 


2 


10 


5.34E+06 


4.79E+06 


5.44EH 


-0(5 


5.12E+06 


o 


20 


1.02E+07 


8.96E+06 


1.05E-I 


-07 


9.15E+06 




15 


1.45E+07 


1.28E+07 


1.52EH 


-07 


1.33E+07 




6 


1.00E+08 


1.08E+08 


7.62E^ 


-07 


7.45E+07 




1 


3.41E+08 


4.40E+08 


1.82EH 


-08 


2.03E+08 


o 


2 


1.99E+09 


2.91E+09 


1.08EH 


-09 


1.21E+09 




1 


1.07E+10 
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we shall need these quantities later. Note that the true He 
yields for set 3 are higher than those of 2a and 2b, in partic- 
ular at the higher metallicities. Note also the lower yields, 
when using mz,=0.01, by a factor of ~2.8. 



4.4 Stellar lifetimes and element timescales 

In the case of set no. 1, 2a and 2b, the stellar lifetimes are 



adopted from Schaller et al. (1992) (Genev a-trac ks) , while 



in set no. 3 they are from Portinari et al. (1998) (Padova- 
tracks). To compare the two sources, lifetimes are shown 
for some selected metallicities and masses in table ^j. It is 
possible directly to compare the stellar lifetimes at solar 
metallicity. Low- and intermediate mass stars seems to have 
shorter lifetimes using the Padova-tracks (between ~40 to 
~70 per cent), while the opposite is true for massive stars 
(but ^10 per cent longer). This age difference may be a 
result of different stellar wind-efficiencies, and the inclusion 
of overshooting in the Padova-tracks. 

By combining our knowledge of stellar yields with the 
respective lifetimes as function of mass, we obtain informa- 
tion on the ejection timescales of different elements. For the 
present purpose, it is the difference between the timescales 
of N and O ejection that is particularly interesting, since the 
time when the ejected N/O has its minimum is equal to the 
time, where the second burst of a pair should appear if max- 
imum scatter is to be obtained. The timescales are found for 
a fixed metallicity at which the ejected masses of O and N of 
one burst are calculated in small timesteps (of the order of 
~ 10 6 yr) until all stars down to 1 Mq have ended their life 
cycle. The time evolution of the element ejection has been 
calculated for set no. 2a and 3 for two metallicities. Fig. [| 
shows the timescales for set 2a. The results for set 3 are very 
similar. 

The production of N in massive stars is different for 
the two metallicities. For solar metallicity, it is released 
much faster for a=l, i.e by more massive stars. This is only 
secondary N though, since the contributors have lifetimes 
shorter than 50 - 80 Myr, hence being more massive than 
8 Mq (compare with the stellar lifetimes given in table |^). 
For set 2a, it is not surprising that the ejection timescale 
of N is much shorter at low metallicity if cv=0.0 (dashed 
lines in fig. ^) than if a=1.0, because no primary nitrogen 
is ejected from intermediate mass stars when cv=0.0. The 
major N production sets in at about 100 Myr at the low 
metallicity. For solar metallicity a large fraction is produced 
within the first 50 Myr, about 40 (2a) to 60 per cent(3). 



100 
ne (Myr) 



Z = 2.0E-i 



100 150 
ne (Myr) 




50 100 150 
Time (Myr) 



Figure 4. The time evolution of O and N ejection after an in- 
stantaneous burst using set 2a, at two different metallicities, and 
two different values of a. Dashed lines correspond to a=0, solid 
lines to a=l. O is to the left and N to the right. The y-axis is the 
mass ejected of N or O, in units of the total mass ejected of N or 
O when all stars have ended their life. The step-like appearance 
of the curves is a consequence of the combined action of discrete 
stellar masses and finite timesteps. t=0 is the time of the burst. 



Almost all of the oxygen is produced in ^30-40 Myr for all 
sets and metallicities. 

Hence, in order to choose the interburst period between 
the bursts of a pair, a minimum value of N/O would appear 
somewhat between 30 and 80 Myr for the low metallicity 
and between 30-40 Myr for the high metallicity. If a choice 
has to be made, 30 Myr would be the best to make, since, 
in general, it is easier to reproduce scatter at low metallic- 
ities due to the relative effect of ejected metals to metals 
already present in the interstellar medium. Following this 
conclusion, the time between the two bursts of a pair is set 
to 30 Myr. 



5 RESULTS 

The model is applied to fit the observations in the N/O-O/H 
plane with the use of the yields in table ^| The accompanying 
results on Y-O/H and O/H vs. gas fraction are displayed as 
well. The order of presentation starts with the results of the 
closed model. 



5.1 The closed model 

The closed model has been applied to all of the yield sets. In 
most plots, the IMF with m_L=0.1 has been used. For those 
employing the 0.01 low-mass limit, it is mentioned separately 
in the text. It is not necessary to do all calculations for both 
IMF's, since the results are the same, except for the gas 
fractions. For all elements, the only difference is a lower yield 
using the low value of mi, hence more massive bursts are 
needed to produce the same abundances as when the higher 
low-mass limit is used. This statement is confirmed in the 
following. 
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Figure 5. Results for a closed model using yield set no. 1 and 
m,L =0.1 Mq The following parameters were adopted: Period 
between pair of bursts is 1 Gyr and the short interburst periods 
are 30 Myr. The mass of the first burst is 1 X 10 6 Mq and the 
masses of the successive bursts are 3 X 10 6 Mq. 
The observations shown are the same as in fig. [l]. 



Figure 6. Results for a closed model using yield set no. 2a and 
mj,=0.1. The following parameters are used: Period between pair 
of bursts is 1 Gyr and the short interburst periods are 30 Myr. 
The mass of the first burst is 1 X 10 6 Mq and the masses of the 
successive bursts are 3 X 10 6 Mq . Finally, the upper mass of hot 
bottom burning, rajjgg, is 5Mq. For this plot, a = 1.1 is used, 
since this value places N/O at the right level. 



5.1.1 N/O-O/H 

Fig. ^ shows the evolutionary path of a double-bursting 
dwarf galaxy, as it has been calculated by the model using 
yield set 1. The point where the evolutionary path starts 
(at 12+log(0/H)^7.1) is where the second burst appears, 1 
Gyr after the first burst. The next point is where the third 
burst would be observed, 30 Myr after the second burst. 
Thus, the upper points of the saw-tooth pattern are cor- 
responding to the first bursts of each pair, and the lower 
points to the second pair-bursts. The scatter is more pro- 
nounced at low metallicities because the mass of O ejected 
relative to the one already present in the ISM is higher, the 
lower the metallicity. It is clear from this figure, why set no. 
1 is used as a reference only. The level of N/O is much too 
high, in particular at low metallicities. The only way one 
can lower the ratio is by assuming a top heav y IM F, since 
the yields from van den Hoek & Groenewegen ( 1997 ) do not 
allow for changing the convection parameters or provide pos- 
sibilities for distinguishing between primary and secondary 
components. Thus, it is not possible to scale the secondary 
components separately to Z=0. Further, neither outflow nor 
inflow is able to cure the problem, in particular not if se- 
lective winds are considered, since this would increase N/O 
even more. 

It is interesting to note that the parameters (convec- 
tion, wind etc. ) use d to fix the yields in van den Hoek & 
Groenewegen (1997) have been found by using a synthetic 
model on AGB stars in the LMC and in the Galactic disc. 
The LMC is not represented in fig. ^, but according to Pagel 
et al. (|1992|) 12+log(0/H)= 8.36 and log(N/0)= -1.22 for 



LMC, which is fitted fairly well after 4 pairs of bursts (but 
may be fitted by 2 pairs of bursts if the mass of each burst 
is larger), thus being consistent. The problem seems to arise 
at lower metallicities. 

Next, set 2a is applied (fig. ^|). The upper mass limit of 
hot bottom burning (HBB) is set to ran bb — 5Mq , compat- 
ible with set no. 3. In this case, a value of a=l.l is necessary 
to obtain the right level of N/O. A larger value of the con- 
vection parameter would produce too much primary N and 
vice versa. Both theoretical predictions and observations 
of AGB stars in the Large Magellanic Cloud are su ggesting 




a level of HBB, roughly co rresponding to a— 2.0 (yan den 



Hoek & Groenewegen 1997, and references therein), which 



Figure 7. Same yield set and parameters as in fig ^j, except that 
mi =0.01, the mass of the first burst is 2x 10 6 Mq and the masses 
of the following bursts are 6 X 10 6 Mq. Compare the results to 
those of fig ^j. It is evident that the observations are explained 
equally well using the low cutoff as when using the higher cutoff, 
if just burst masses are about twice as massive. 



is well above the value obtained here. If a really is so large 
in AGB stars in the LMC, an explanation might be that the 
value of the convection parameter depends on metallicity, 
since the metallicity of the LMC is high (about half-solar). 
A higher value of a at high metallicities is not ruled out by 
the model as long as it is equal to 1.1 during the first few 
bursts. 

A large spread is evident and sufficient to explain the 
observed scatter. Here, the plusses, representing the DLA- 
systems, are neglected. A similar run is presented in fig. ^, 
only employing the IMF having m_L=0.01 Mq. As expected, 
it is seen that the result resembles the situation in fig. |^, if 
one assumes all bursts to involve twice the mass of the bursts 
in fig. H 

Figures ^ and ^ give a good comprehension of the out- 
come of the model, but some considerations should be made 
before any solid conclusions are drawn. Firstly, abundances 
are measured in H n-regions only, and therefore only at the 
'positions' of the bursts, not in between. Secondly, the evo- 
lutionary path represents the evolution of one dwarf galaxy. 
Another dwarf galaxy, having different parameters such as 
total mass, burst masses etc., would have another evolution, 
and thus another path in the N/O-O/H plot. Thus, fig. g 
shows the evolutionary paths of dwarf galaxies having differ- 
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Figure 8. Three closed models, each with different burst masses, 
using set no. 2a. All other parameters are the same as in fig. fil 
For each burst mass, two lines have been drawn: one through the 
upper points of the saw-tooth pattern and one through the lower 
points. The dashed lines correspond to burst masses equal to f 
per cent of the total mass of the galaxy, the solid lines to 3 per 
cent and the dotted lines to 5 per cent 



ent fractional burst masses. Yield set 2a has been used, and 
all of the parameters, except the mass of each burst, are the 
same. To make the plot more clear, lines have been drawn 
through the first bursts of each pair, and through the sec- 
ond bursts. The former lines are almost coincident, because 
the time interval between two pairs of bursts is ~ 1 Gyr, 
enough to release almost all of the nitrogen. Thus, larger 
bursts increase both N and O, adjusting N/O to be roughly 
unchanged. The case is different for the lines through the 
second pair-bursts, because the N from the first pair-burst 
has not yet been released after 30 Myr. Thus, a larger burst 
mass gives an increase in O abundance only. 

Two results may be deduced from fig. ^ Firstly, the 
scatter of the observations is perfectly explained, because 
an observational scatter of ~0.1 dex has to be taken into 
account. Furthermore, the area between the lower line cor- 
responding to 5 per cent burst masses and the upper line 
is 'filled' in the sense that the distribution of burst masses 
may be everything between zero and 5 per cent to match 
the observations. Thus, the distribution of abundances is 
accounted for by using the models of dwarf galaxies having 
various reasonable relative burst masses. The other result is 
the slight increase of the evolutionary path toward higher 
metallicities. However, the small upturn may not be a con- 
sequence of dominating secondary N production, but rather 
the decreasing O yield at higher metallicities, see table |^. 

Set no. 2b is different from no. 2a only in that a modest 
mass-loss from massive stars is assumed. The effect of this 
change on the evolutionary path is seen by comparing fig. 

with fig. ^. It is possible directly to compare the model 
outputs of the two figures, since all parameters, including 
the burst masses, are identical. When using set 2b instead 
of 2a, the N/O ratio has an even smaller upturn at higher 
metallicities, because the O-yield is higher. Thus, even when 
N production is increasing at higher metallicities, so is the 
oxygen, keeping the ratio down and thereby hiding the sec- 
ondary behavior of N. Still, the conclusions are the same: 



Figure 9. Evolutionary paths of three closed models using yield 
set 2b. The lines are drawn as in fig. ^ with identical parameters 
so the dashed lines correspond to burst masses equal to 1 per cent 
of the total mass of the galaxy, the solid lines to 3 per cent and 
the dotted lines to 5 per cent 




8.0 

12 + Log(0/H) 



Figure 10. Three closed models each with different relative burst 
masses using the yields of set 3 (Padova). The meaning of the lines 
is the same as in fig. |i[ The dashed line corresponds to burst 
masses equal to f per cent of the total mass of the dwarf galaxy, 
the solid line to 3 per cent and the dashed-dotted line to 5 per 
cent. It is worth noting that the mass interval for producing pri- 
mary N is the same as in fig. ^. For these runs, all primary yields 
have been extrapolated for metallicities lower than the range cov- 
ered by the yield sources. 



both the level of N/O and the scatter is explained, indepen- 
dent of m,L ■ 

Using yield set no. 3, the same kind of plot is con- 
structed as those above. Here there is no explicit assump- 
tion on the convection parameter, but for reasons explained 
below, the primary component of N (and O) for Z<0.008 
has been extrapolated linearly, using the yields of Z=0.008 
and Z=0.02 instead of just adopting the primary yield at 
Z=0.008. The run of the model is presented in fig. |Io[ At 
least two features have to be mentioned. The first is that the 
scatter is somewhat smaller than for set 2a (or 2b) and the 
second is an upturn in N/O at higher metallicities, being 
more pronounced, than it was for set 2a and 2b. Unfortu- 
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Figure 11. An evolutionary track, using set 3 yields. All burst 
masses are 3 X 10 6 Mq. The yields of primary components for 
metallicities lower than Z=0.008 are set equal to those given for 
Z=0.008. 



nately, the upturn is taking place for higher oxygen abun- 
dances than those of the observations. 

The primary component in set 3 has also been mod- 
elled in another way: For the same choice of parameters 
as in fig. [l^ (solid line), the model has been applied with 
primary yields equal to the one at the lowest metallicity for 
which the yields are given (Z=0.008), with the resulting evo- 
lutionary path plotted in fig. [n]. A comparison shows that 
an adoption of the primary yields at Z=0.008 will produce 
too little primary N. It seems to be necessary to make pri- 
mary production more effective, the lower the metallicity, 
as it is done by the linear extrapolation, because the pri- 
mary N yield at Z=0.008 is slightly larger than it is for solar 
metallicity. This is actually an important result, because it 
provides a constraint on the produced primary N, indepen- 
dent of stellar models. The explanation for this result might 
be that the number of thermal pulses in the AGB phase of 
low metallicity stars is higher, thus converting more C and 
O into N via the CNO-cycle. Hence, in the following, all 
model runs using set 3, will use linear extrapolations of pri- 
mary N yields, below Z=0.008. Of course, one could make 
the IMF steeper thereby producing less oxygen. However, 
observations do not support a steeper IMF. 

When comparing the outputs of the model, using set 
2a (2b) and 3, it is important to remember that in set 2a 
and 2b, the lowest metallicity for which yields are available 
is Z=0.004, at least in the critical, and very important mass 
interval 4 to 5 Mq, while it is as high as Z=0.008 in the same 
mass interval for set 3. Thus, whether primary N yields at Z 
<0.008 are extrapolated from Z=0.02 over Z=0.008 or are 
adopted to be equal to the yield at Z=0.008, the result is 
quite uncertain, and stellar yields at much lower metallicities 
are definitely desirable. 



5.1.2 Y-O/H 

It is worth emphasizing that the helium abundances are fit- 
ted with input parameters not different from those used to 
fit the N/O ratio. For the closed model presented in fig. ^ 
by the solid line, Y has been calculated as well. Y against 



Figure 12. The evolution of the helium mass fraction Y as a 
function of O/H, using yield set no. 2a and mj,=0.1, plotted as 
a solid line. The model output of Y and O/H are from the same 
run as marked by the solid line in fig. || thus having a burst mass 
of 3 per cent of the total galaxy mass. The meaning and sources 
of the symbols are the same as in the N/O - O/H plots, except 
that the DLA systems are not represented. Note that the O/H 
axis is not logarithmic. 



O/H is shown in fig. The evolutionary Y path does not 
show any scatter, though the dwarf galaxy is double burst- 
ing, because about half of the helium is produced in massive 
stars, together with the oxygen. The fit is very good and 
well within the error bars. For the two most metal deficient 
objects in the plot, IZwl8 and SBS0335-052, it seems to be 
necessary to assume smaller burst masses, say 1 per cent of 
the total mass - the model with the dashed lines in fig. ^[ 
Thus, these objects can be explained by experiencing their 
second burst and having burst masses equal to 1 per cent of 
their total mass. The output of the model gives the abun- 
dances just before every burst. Hence, to calculate the model 
value of dY/dZ, the Y,0/H values were fitted using a lin- 
ear least-square fit, assuming Z=20(O/H). The result is 1.0, 
lower than the (very uncertain) value of 2.6 derived from the 
observations, see section The corresponding result using 
set 2b is not shown, since it resembles the one by set 2a. Cer- 
tainly, it is true that the He yield is smaller for modest mass 
l oss, but only for very massive stars and high metallicities 
( |Maeder 199% . 

The dependency of He on O is more pronounced, when 
using the Padova yields, as seen in fig. [li| Calculating 
dY/dZ, using this yield set, in the same way as done for 
set 2a one obtains 3.1, higher than the value calculated for 
the observations (2.6). From table ^ it is clear that the true 
He yields for set 3 are a factor of ~3 higher than those of 
sets 2a and 2b, but it is unclear why the Padova tracks pro- 
duce that much helium. This feature makes the fit not quite 
as good as it was for set 2a (2b). However, the p rimordial 
helium determination is only certain within 0.003 (Izotov et 



al. 1997a). Setting the primordial He abundance equal to 



the lower limit (0.240) in the initial conditions of the model 
does not change the N/O-O/H evolutionary path, but it 
does lower the level of Y in fig. [ll]. The new evolutionary 
path of Y has the same slope as the solid one of course, but 
sufficiently lower to be a fairly good fit, as it appears well 
within the observational uncertainties. 
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Figure 13. The evolution of Y using the Padova yields. The solid 
line is the evolutionary path, using a primordial He abundance 
Yp=0.243, from the same model run as used for the solid line in 
fig. [w| and the dashed-dotted line is the corresponding path, as- 
suming Yp=0.240, still within the uncertainties of the primordial 
He determination. The mass of each burst is 3 per cent of the 
total mass in both cases. 



5.1.3 O/H-n 

Gas fractions may be useful in serving as a further con- 
straint on chemical evolution models. It has been pointed 
out by several authors that closed models, such as the one 
presen ted above, are unable to explain the observed gas f rac- 
tions (|Mattcucci fc Chiosi 1983| ; |Matteucci fc Tosi 1985] |Ca 



rigi et al. 1995). Thus, gas fractions have been calculated 
just before every burst using eqs. [23] and ^ allowing for 
comparison with some observed values from the literature. 

Problem lies in the interpretation of the observations. 
The total dynamical mass of a dwarf galaxy also includes 
dark matter (DM). This dark is unlike to participate in the 
chemical evolution. Our model is including baryonic mat- 
ter only. Thus, in order to compare the observations with 
the calculated gas fractions, one has to be sure that no sig- 
nificant amount of non-baryonic DM is included in the to- 
tal ma ss estimates. However , as pointed out by s e veral au- 
thors (|Bi-inks & Klein 1988[ fcumai fc Tosa 1992| |Carigi et 



al. 199£; MeureT*1998j ) dwarf galaxies, including BCGs, are 



DM-dominated. Consequently, dynamical estimates of total 
masse s are u seless, when considering dwarf galaxies. Carigi 
et al. (1999) tried to solve the problem by using dynamical 



estimates within the visible Holmberg radius, to avoid in- 
clusion of a DM halo. Unfortunately, this leaves us with a 
very small sample of observations, and still the amount of 
DM within the Holmberg radius is uncertain. 

However, gas fractions are one of the few possible ways 
of constraining chemical evolution models, so a discussion 
of gas fractions, using the same models as in the preceding 
paragraphs should be included. To avoid erroneous conclu- 
sions, we used lower and upper limits of the gas fractions. 

The gas fractions used are presented in table ^. The 
gas masses included are Hi-masses multiplied by 1.3 to 
account for the helium content. The amount of H2 is ig- 
nored. All dynamically estimated gas fractions have been 
scaled to a Hubble constant Hq = 65kms _1 Mpc _1 using 



M gas oc (distance) 2 and Mtot oc (distance), so the gas frac- 
tion fi oc (distance) oc Hq . 

DM is implicitly included, and molecular hydrogen is 
ignored, hence these gas fractions must be regarded as 
lower limits. Unfortunately, even the lower limits seem to be 
very uncertain. For instance, the gas fraction of IZwl8 was 
first calculated using the data from Staveley-Smith, Davies 
& Kinman (1992). However, after He-correction and Hq- 
scaling, the gas fraction became larger than 1. This is also 
the case for some of the galaxies in the Thuan & Martin 



( 1981 ) sample (not included in the present sample). In the 
case of IZwl8, the value 15 pe r cent gas has been adopted 
from Matteucci & Chiosi (1983). Extreme upper limits have 
been obtained assuming that each of the galaxies are ex- 
periencing their first burst, and are observed during max- 
imum luminosity. Indeed, this is an extreme upper limit, 
and should be regarded as such. Adoptin g the evolutionary 
stellar models of Leitherer & Heckman (1995), the M/Lb 
ratio may be estimated for a burst, M being the mass of 
stars. Leitherer & Heckman calculated the blue luminosity 
for a 10 6 A/q starburst as a function of time. To be consis- 
tent, their Salpeter IMF has been scaled to ours. Indeed, 
two different IMF's have been used so far, but to be sure 
that upper gas fraction limits are obtained, the high lower 
limit, r«L=0.1Mg, is used, as this gives the lowest M/Lb, 
hence the highest gas fraction. The scaling is very simple, 
since the only difference is to scale the lower mass limit 
from 1 to 0.1 Mq. Hence, all luminosities from Leitherer & 
Heckman are multiplied by 0.39. From their fig. 9 (Z=0.1 
Zq), the minimum magnitude is read off to be -16.4 corre- 
sponding to a maximum luminosity ~ 5.6 x 10 s Lq using 
the absolute blue solar magnitude 5.48. After IMF scaling 
this gives ~ 2.2 x 10 8 i Q , so M/L B ~ 0.004. In table | 
the observed blue luminosity, scaled to h=0.65, is given for 
each galaxy. Table ^ gives the corresponding gas fraction, as 
fi = Mgas/Mgaa + M s t a rs- If IZwl8 is experiencing its first 
burst, this upper limit should be realistic as its true gas frac- 
tion. The large difference from the dynamically deriven gas 
fraction is caused by a combination of influence from DM, 
poor total dynamical mass estimation or/and that IZwl8 is 
not experiencing its first burst. 

To constrain the models more, three more galaxies have 
been included, for which we know that they are definitely 
not experiencing their first burst. One of them is the LMC. 
Strictly speaking, the LMC is not a dwarf galaxy, but in 
many ways it behaves as such. For these three galaxies, 
M/Lb is assumed to be 0.5, not as extreme as above, thus 
obtaining upper limit gas fractions lower than for the other 
galaxies. Does this assumption still provide an upper limit 
for the gas fractions? The answer is confirmatory, if the as- 
sumed M/Lb is a lower limit. The rati o is estimated us- 



ing the star for mation histories of LMC ( Geha et al. 199£ , 
fig. 7(c)), SMC (|Pagel fc Tautvaisene 1998| ), and NGC6822 
( |Gallart et al. 1996 , fig. 12). The general procedure is quite 
simple. First we find Lb of the objects by using a stellar 
population synthesis model and t he SFH. The blue magni- 
tudes in Chariot & Bruzual ( 1991, fig. 5, dash ed line) are cal- 
ibrated by using Leitherer & Heckman (1995, fig. 10, 0.25Zq, 
solid line). When we know the age of a starformation epoch, 
we may find the blue magnitude, hence the luminosity by 
reading off figure 5 of Chariot & Bruzual. The results are 
valid for m,£ = lM0, mu~100MQ and a starformation rate 
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Table 5. The data used for gas fraction plots. 



Name 


Lb (Lq) 


M gas (10 8 Mq) 




logfi 


12 + log(0/H) 




(1) 


(2) 


(3) 


(4) 


(5) 


Mrk 5 


0.35E+9 


2.25 


-2.65E-3 


-0.13 


8.21 


36 


0.96E+8 


0.33 


-5.03E-3 


-1.19 


7.85 


67 


0.17E+9 


0.28 


-1.07E-2 


-0.52 


8.21 


178 


0.96E+8 


0.24 


-6.82E-3 


-0.09 


7.95 


71 


0.11E+10 


17.1 


-1.16E-3 


-0.01 


7.89 


209 


0.76E+8 


0.87 


-1.52E-3 


-0.38 


7.81 


600 


0.44E+9 


5.71 


-1.34E-3 


-0.68 


8.01 


IZwl8 


0.10E+9 


1.12 


-1.63E-3 


-0.83 


7.22 


IIZw40 


0.23E+9 


7.79 


-5.04E-4 


-0.29 


8.14 


IIZw70 


0.15E+10 


5.88 


-4.31E-3 


-0.54 


8.06 


VIIZw403 


0.53E+8 


0.66 


-1.40E-3 


-0.39 


7.73 


IZw49 


0.17E+10 


7.96 


-3.76E-3 


-0.28 


8.03 


IZwl23 


0.36E+9 


0.85 


-7.31E-3 


? 


8.10 


CG 1116+51 


0.15E+9 


3.63 


-7.00E-4 


-0.07 


7.52 


SBS 335-52 


0.98E+9 e 


1.94 e 


-8.80E-3 


? 


7.33 


NGC 6822 


0.18E+9 a 


2.53 a 


-0.13 d 


-0.80 


8.19 


SMC 


0.99E+9 6 


4.66 c 


-0.31 d 


-0.38 a 


8.04 


LMC 


0.47E+10 6 


6.90 c 


-0.65 d 


-0.92 a 


8.36 



(1) The blue luminosities are from Thuan & Martin (1981) with four exceptions, see (a,b,e). 
The luminosities are scaled to h=0.65. (2) The gas masses are from Thuan & Martin, though 
see (a,c,e). All values are scaled to h=0.65. (3) The logarithm of the calculated upper limit gas 



to do their very first burst, see text and note (d). (4) The logarithm of the observed, lower 
limit gas fraction, fi is the gas mass divided b y the dynamically estim ated t otal mass from 
Thuan & Martin (1981), Matteucci & Chiosi (1983) and Taylor et al. (1994). The question 
marks indicate that the gas fraction is not found, so in these cases, only the upper limits are 
used. All values are scale d to h =0.65. thoug h see note (a). (5) The o xygen abundances from 
Kobulnicky & Skillman ( 1996 ), Pa gcl et al. ( 1992 ) and Izotov et al. ( 1997a ). 

(a) Values are from Lequeux et al. (1979), obtained by h-independent distance measurements. 
Hence, these values are adopted without any scaling. 

(b) These luminosities are calculated from the absolute magnitudes given in Matteucci & 
Tosi (1985), using the solar absolu te blue magnitude 5.48. 

(c) Values from Carigi et al. (199E ). They give the gas masses directly for h=0.65 and account 
for i?2- To be consistent with the rest of our sample, their values are divided with 1.1 to 
obtain gas masses not including molecular hydrogen. 

(d) These galaxies are known for sure not to experience maximum luminosity of their first 
starburst at the present epoch. Thus, more realistic upper gas fraction limits have been 
calculated, assuming their M/Lb ratio to be 0.5, see text for further details. 

(e) Values from Thuan, Izotov & Lipovetsky (1996) 



of IMoyr -1 , since the magnitudes are calibrated by using 
Leitherer & Heckman. By doing this we ignore the very few 
stars with masses between IOOM0 and 125M© originally 
included by Chariot & Bruzual. This will not have any sig- 
nificant effect on the results. 

For LMC we have the SFR in relative units. We have 
splitted the SFH into two star formation epochs, one that 
started 2 Gyr ago and one 12 Gyr ago. Integration of the 
SFR and weighting the two read-off luminosities with re- 
spect to SFR leads to M/Lb =0.4. By scaling the luminosi- 
ties to ttil=O.1M0 and m_L=O.OlM0, we obtain Af/Ls=1.0 
and 2.4, respectively. 

For NGC 6822 we are given the SFR in absolute units, 
but we restrict ourselves to use the relative SFR only. Again 
we split the SFH into two epochs. A very recent one that 
started 200 Myr ago and a very old one that started 15 Gyr 
ago. As we ignore the epoch that stopped 5 Gyr ago, we 
will obtain a lower estimate of M/Lb- Using the method 



outlined above, we get M/Ls=0.16, and 0.42, 0.96 for the 
0.1 and 0.01 lower mass limits, respectively. 

The lack of knowledge of the SFH in the SMC is re- 
markable. Hence, we h ave u sed the SFH from the model by 
Pagel & Tautvaisene 



199c 



which has succes in explain- 
ing e.g. [Fe/H]. The SFR is given in absolute units in their 
table 2, so integration yields the total mass of stars. Using 
L B = 0.99 x 10 9 Mq from table |, we arrive at M/L s =0.43. 

In conclusion we have that M/Lb =0.5 is a conservative 
lower limit, hence giving an upper limit on the gas fraction. 

The results from the calculations are shown for yield 
set 2a (top plots) and 3 (bottom plots) in fig. [l4|. In both 
cases, using the IMF with rriL = 0.1 Mq, the fits are al- 
ways between the lower and upper limits, except for the 
three 'moderate upper limit' galaxies. Lowering the yield by 
adopting rriL = 0.01 Mq, helps somewhat, but only one of 
the three moderate galaxies is fitted, namely NGC 6822, the 
two others being out of range. 

Since the model fails in explaining the more restrictively 
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Figure 14. 12+log(0/H) as a function of gas fraction, fj, = 

■jjj — M +m t ' ^ P plots are for set 2a, for mi = 0.1 Mq to 

the left and mi = 0.01 Mq to the right. Bottom plots are for set 
3, again mi =0.1 Mq to the left and Mi ow = 0.01 Mq to the 
right. The dashed lines are the results of our closed models, using 
burst masses equal to 2 X 10 6 Mq for the plots on the left and 
4 X 10 6 Mq for the plots on the right. The stars are the lower 
and upper limit gas fractions for the majority of the sample, con- 
nected by lines. The diamonds are the limits given for the three 
galaxies, having more moderate upper limits. 



chosen galaxies, the yields have to be lowered further. The 
number of low-mass stars could be increased, but it would 
be parameter-gambling, not appealing very much to a phys- 
ical control of the model. However, it should be noted that 
Carigi et al.( 1999 ) use this possibility to solve the problem. 
Keeping the IMF, the only way of lowering the yield is to 
open the model and allow for galactic win ds, ordinary or 
enriche d (Sarigi et al. 199E| ) or/and inflow (Pagel & Taut- 
vaisene 1998). 



5.2 The model including enriched winds 

From this place forth, the model is opened, allowing for gas 
exchanges with the intergalactic medium. All simulations 
have used mi = 0.01 Mq. 

The results of including enriched winds are shown in 
fig. ^| for set 2a and fig. [b| for set 3. As seen, the wind 
efficiencies employed are the lowest acceptable for fitting 
the gas fractions of the Magellanic clouds. Using the same 
model, the outputs have been plotted with the observations 
for both N/O-O/H and Y-O/H. 

However, removing that much oxygen raises the N/O 
ratio high above the level of the observations. One could try 
to reduce N/O by lowering the value of a, in the case of set 
2a, but a value much lower than 1, seems to be quite unre- 
alistic, seen in the light of recent work on stellar evolution, 
e.g. Marigo fll998p . 

For set 3 (fig. it is clearly seen that setting the 
primordial He fraction equal to its lower uncertainty value, 
as done when using set 3 in the closed model (see fig. |l^) , is 
not capable of explaining the He fractions, when introducing 
these high- efficient enriched winds - the fits are much too 
poor. 

Hence, introducing enriched winds makes it impossible 




:.kl 




Figure 15. The results of models with enriched winds. The plots 
are made, using an efficiency factor equal to 0.8, adequately to 
fit the gas fractions. The yield set in use is 2a. The values of 
Mhbb and a are the same as used for the closed models, thus 
being 5 Mq and 1.1, respectively. All burst masses are equal to 
6 X 10 6 Mq. 



to fit gas fractions and N/O simultaneously. Further, in the 
case of set 3, problems also arise in fitting the He fractions. 



5.3 The model including ordinary winds 

Though keeping the mass of each burst as a free parameter, 
it turns out that 2 x 10 6 Mq is fulfilling the requirement on 
placing the second burst close to the abundances of IZwl8. 
Thus, it is appropriate to keep this burst mass constant in 
the following discussion. The results are shown, in fig. [l7|for 
set 2a and fig. |l^ for set 3. 

It was shown above that the scatter in N/O could be 
explained by our model. Hence, the scatter will not be given 
much attention in the following, as the work will be concen- 
trated on getting the right level of N/O, coincident with the 
gas fraction fitting, if possible. 

From the top plots, it is evident that a mass removal 
of 5 times the burst mass does not appreciably change the 
appearance of the fits from the closed model. Actually, the 
O-yield is lowered only by ~0.3 dex as seen when comparing 
to fig. JO], using 2a yields. The point may be that we are deal- 
ing with a discontinuous star formation. Note for instance 
the upper linear part of the model curve. The increase in 
O/H along this part is due to one burst only, as seen when 
comparing to the corresponding N/O plots. The explanation 
is that the smaller the absolute gas mass, the larger the ef- 
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Figure 16. The results using the same parameters as in fig. |15| 
but with yield set no. 3 (Padova yields). The efficiency factor is 
0.7, chosen to fit the gas fractions. 



feet of oxygen enrichment. Thus, even when adopting higher 
wind efficiencies, it is not possible to bring the oxygen yield 
sufficiently down. The same mechanism is responsible for 
increasing the scatter in N/O for higher wind efficiencies. 
Hence, it is possible to produce large scatter without having 
a large burst mass. 

Note also that the evolutionary tracks stay close to a 
gas fraction of 1 for a relatively large part of the evolution. 
A simple calculation will show, how difficult it is to bring 
the gas fraction down. For instance, assume the first burst 
to turn 2 x 10 6 Mq into stars, and set Wism = 5. Then, just 
before the second burst, the gas fraction will be 



/'• : 



Mgas 

Mtot 



M gas (0) — 2 x 10 6 Mq — W 7 M e 



M tot (P) - W 7 M e 



(30) 



0.' 



since the initial mass of the dwarf galaxy is M ga s(0) = 
Mtot(0) = fO s M0. Using the equation again, one finds 
/i = 0.95 just before the third burst. These values are lower 
limits (!), since the increase in gas from stellar ejecta is dis- 
regarded. 

Both the level of N/O and slope of Y are satisfied. Note 
the similarity between the Y plots presented here and those 
presented by the closed model. This is not surprising, as 
about half of the helium is produced by massive stars, hence 
following oxygen. Status is that all observations are matched 
reasonably well by including ordinary winds, but no better 
than for the closed model. In particular, it is found that it 



Figure 17. The results of using a model including ordinary 
winds. The plots are made using Wism = 5. Note the increasing 
scatter in N/O as a function of Wism explained in the text. The 
yield set in use is 2a. The values of Mh b b and a are still 5 Mq 
and 1.1, respectively, and the mass of each burst is 2 X 1O 6 M0. 



is impossible to bring the yields sufficiently down to explain 
the gas fraction intervals of the Magellanic clouds. 



5.4 Comparing the ordinary wind model with an 
analytical model 

A few calculations have been performed, using a simple an- 
alytical model, employing continuous star formation and 
Hartwick-outflow. The outputs of the analytical model and 
the numerical model are then compared. The comparison is 
only performed for yield set 2a, since the inclusion of winds 
is identical using set no. 3. 

The general formulation of Hartwick-outflow is 



dM d v 



-Wism 



dM s 



(31) 



dt dt 
where Mdw is the mass of the dwarf galaxy and M s is the 
mass of stars. The next step is to use 



dZ 
~dM~ s 



M 



(32) 



obtained by considering the changes in Z as a result o f stel- 
lar ejecta, star formation and outflow, see e.g. Pagel ( 1997 , 
his eq. 7.37). Instantaneous recycling is assumed, and so is 
the absence of inflow. Separating the variables and using 



5M„ 



-8M S - Wism8M s , one finds 
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Figure 18. The results using the same wind and mass parameters 
as in fig. nil but with yield set no. 3. 




-0.5 
Log y. 



Figure 19. The comparison between the numerical bursting 
model employing ordinary winds with Wj$m = 5 and burst masses 
2 X 10 6 Mq, and the continuous star forming analytical model. 
The solid line is the numerical bursting model. 



1 + Wism 



■In 



M g (Z = 0) 
M a 



(33) 



where the yield is assumed constant. The total mass of the 
dwarf galaxy is written 



M dw = M d w{Z = 0)-WismM s 
= M g {Z = 0)-WisM(M d « 



M 



(34) 
(35) 




Figure 20. Comparing the analytical model from fig. |h| with the 
semi-continuous numerical model. The solid line is the numerical 
model and the dashed line is the analytical model. See text for 
details. 



Isolating M g (Z = 0) and inserting it into eq. p3| gives 
1 + Wism 



exp (f (1 + Wism)) + Wism 



(36) 



exp 



1 + Wism 



V ) exp (Wismj) + Wism exp (-|) 



where fj, is the gas fraction 



It is seen that for a closed 



model, (i = exp ( — ~J. For our numerical model, the yield 

is metallicity-dependent. Hence, exp (— , equal to the gas 
fraction at Z for a characteristic value of p, is found from 
the closed numerical model. Inserting this into eq. |3(| finally 
gives the gas fraction of the analytical model, now employing 
metallicity-dependent yields and outflow. 

The results are shown in fig. [is] for Wism=5. The nu- 
merical fit (solid line) is the same as in fig. [n]. Note that for 
a large part of the evolution, the two models are coincident, 
but at the last part they are differing from each other. It is 
likely that the difference is a result of using instantaneous 
bursts instead of a continuous SFR. The point may be that 
the mass of gas is small at the last part of the evolution. 
Thus, the relative amount of ejected oxygen from one burst, 
is high compared to the gas mass, hence giving a large in- 
crease in O/H. 

To make a further check on this hypothesis, a contin- 
uous version of the numerical model has been made. This 
is done in an approximate way, introducing the following 
changes to the above numerical outflow model: 

(i) The bursts are all single bursts, i.e. all interburst pe- 
riods are the same. 

(ii) The interburst periods are made short. The shorter 
the interburst period, the more continuous the star forma- 
tion. An interburst period of 3 Myr was sufficient for our 
purpose. 

(iii) The masses of the bursts are calculated by assuming 
a mean SFR. The calculations here assume < SFR >=0.002 
MQj/r -1 , corresponding to one burst of mass 2 x 10 6 Mq 
every Gyr. A reasonable value compared to the calculations 
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above. The masses of the bursts are then equal to 0.002 
M Q yr' 1 - 3 Myr = 6000 Mq. 

All other calculations are performed exactly the same 
way, as they were, using the double bursting model. The 
calculations are performed for Wism~ 5 and compared to 
the analytical model from fig. From fig. it is seen that 
the two models are almost identical. The small offset may be 
caused by the calculation of the analytical model, adopting 
the yields from the closed numerical bursting model. The 
similarity between the two models confirms the suspicion 
that instantaneous bursting models behave differently from 
continuous models, as the numerical model is the same that 
produced the fit in fig. [lS^, except for the semi-continuous 
star formation. 



5.5 The model including ordinary winds and 
inflow 

The introduction of inflow gives two new parameters: the 
total mass accreted onto the galaxy Mo and the duration 
of the inflow event Tin/. Only results for a few parameter 
choices are shown and interpreted by close examination of 
the results. 

The gas fractions of the three selected galaxies provide 
strong restrictions on the possible parameter space, since 
inflow and outflow balance each other for obtaining the right 
gas fraction values. Hence, if one increases the rate of inflow, 
one has to increase the rate of outflow as well. Otherwise, 
the modelled gas fractions would be too high, and the three 
selected gas fraction intervals would never be reached. 

It is assumed that the mass of the original gas cloud is 
at t=0. However, inflow of gas increases the mass of the 
gas cloud, until it starts to form stars at a mass threshold 
M(t'). The time t' is calculated using eq. EQ, giving 



Ml/ 



In 



Mo 



Mo - M(t') 



(37) 



For consistency with previous calculations, the threshold is 
assumed to be 10 8 Mq. Thus, after a time t' , the first burst 
appears. As for all the previous models, the first burst is 
a single one, and all other bursts appear in pairs. Hence, 
the appearance of the first burst is delayed by the time t' , 
when comparing to the models not including inflow. This 
time difference is included, when calculating the age of the 
system. Remember that the model terminates if the system 
is older than 15 Gyr. 

For yield set 2a, the results are displayed in figs. |2l], ^ 
and figs. 

For the parameter values chosen very small scatter re- 
sults in the N/O diagram in fig. The plots were made us- 
ing a high inflow/outflow rate to lower the yield sufficiently 
to fit the gas fractions of NGC 6822 and the Magellanic 
clouds. The price to pay seems to be a very small scatter 
in N/O. Hence, in fig. |2^ a low rate of inflow has been as- 
sumed, in particular at late stages of evolution, as Ti n f is 
decreased to 2 Gyr. Now, it is possible to explain the scat- 
ter in N/O, but the result is an unsatisfactory gas fraction 
fit. Obviously, a large inflow rate of primordial gas buffers 
the impact of O-rich ejecta on the ISM abundances, result- 
ing in small scatter. 

To investigate further, whether it is possible or not to 




Figure 21. The results of including both inflow and ordinary 
winds. Here the yield set used is no. 2a. The plots show a run 
using Ti„f=5 Gyr and Mo = 8 X 10 8 Mq. The burst masses are 
3 x 10 6 M e , and W IS m=8. 



explain the scatter in the N/O-O/H diagram and gas frac- 
tions simultaneously, the outflow parameter is kept at 8, but 
the burst masses are increased to 8 x 10 6 Mq to give a larger 
scatter. To balance the high outflow rate (remember that the 
outflow rate is proportional to the burst mass) a rather high 
inflow rate is needed, (ri„/ ,Mo) — (5 Gyr,10 9 Mq). 

The results are shown in figs. and Note that the 
high outflow/inflow rate forces N/O to increase dramatically 
with increasing O/H, because inflow of primordial gas de- 
creases O/H, but not N/O. At the same time, the strong 
outflow decreases the O abundance more than the N abun- 
dance, since the composition of the ISM just before each sec- 
ond pair-burst is O enhanced, hence giving an O enhanced 
wind (though defined and treated as an ordinary wind). 

The corresponding gas fraction plot fits NGC 6822 and 
SMC, but the LMC gas fraction and abundance are never 
reached. Note the 'sawtooth' behavior in the gas fraction 
plot. During the long interburst period, the high inflow rate 
deposits a large amount of H in the ISM, but no O, hence 
decreasing O/H until the next burst enriches the ISM again. 
However, the gas fraction does not increase, because the in- 
flow is balanced by the high outflow rate. From these last 
results, it seems to be possible to explain both the scat- 
ter in N/O and the gas fractions simultaneously, but note 
that in our struggle to explain the gas fractions of the three 
selected galaxies, some 3-4 other dwarf galaxies are not fit- 
ted within the intervals. Note also that the closed model 
explained those gas fractions quite well! 
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Figure 22. Still using set 2a, though trying to increase the degree 
of scatter. Now, the two parameters of inflow are r in f= 2 Gyr and 
M = 2 X 10 8 M Q . The burst masses and wind parameter have 
the same values as in fig. EU 




7.0 - 
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Figure 23. Another model using yield set 2a, employing burst 
masses equal to 8 X 10 8 M Q , Mo = 10 9 M Q , r inf =5 Gyr and 
Wism = &- Note the 'sawtooth' behavior caused by the high inflow 
rate. 



The results using yield set 3 are displayed in fig. |2^. The 
gas fractions of the three well-known galaxies are properly 
fitted, though the same remark as above has to be made, 
namely that the fits are outside the intervals of 3-5 other 
dwarf galaxies. 

The right level of N/O is obtained. Note in particular 



Figure 25. Results using yield set 3. For the model represented 
here the adopted burst masses are 6 X 10 6 Mq and M} = 7 X 
10 8 M Q . Finally, W IS m=5 and T inf =5 Gyr. 

the increasing N/O in the left plot, following the trend of 
the observations. As for set 2a the scatter in N/O is rather 
small, and insufficient to explain the observed scatter. 

The four Y-O/H plots have not yet been discussed, but 
it may be done in a few words. The model predictions match 
the observations quite well, but they are of course still show- 
ing the difference in slope between set 2a (or 2b) and 3. 
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6 DISCUSSION 

Numerical models calculating the chemical evolution of gas- 
rich dwarf galaxies have been presented. The models have 
been fitted to a sample of abundance and gas fraction ob- 
servations. A chemical evolutionary model has to fit these 
observations simultaneously. It is clear that a realistic model 
is not just a couple of equations including a lot of parameters 
that one can change until all data are fitted. One has to re- 
member that the equations are applied to a physical system, 
obeying physical laws. Hence, before doing any calculations, 
some considerations were made about ejecta dispersal and 
mixing processes. It was found that the processes involved 
are complicated and no complete theory exists yet. How- 
ever, the four most important features in the enrichment 
process are found to be stellar ejecta, giant Hn-regions, 
wind-driven superbubbles and SN-driven supershells. Wind- 
driven superbubbles arise shortly after a burst due to the 
strong stellar winds of massive stars. The radii of these su- 
perbubbles are in general smaller than the radii of the H II- 
regions. As the massive stars explode as SNe, a supershell 
is swept-up, soon catching up on the superbubble. However, 
a simple calculation showed that their radii become compa- 
rable to the radii of H n-regions only in the late stages of 
Hn-region existence. Hence, it is suggested that observed 
emission lines, used for abundance determinations, arise in 
the ionized medium outside of the superbubble/supershell. 
Numerical hydrodynamical models indicate that SN-ejecta 
always stay within the supershell. If this is true, the observed 
abundances are not affected by the SN-ejecta from the stars 
producing the Hn-region, hence being typical for the ISM. 

The details of supershell evolution are still not known, 
but both theory and observations support star formation in 
expanding shells. From theory, it is expected to start not 
earlier than about 20 Myr after the starburst, resulting in 
star formation involving a mass comparable to the mass of 
the burst that initiated the supershell, hence appearing as 
a double burst. The double bursting mode of our numerical 
model ensures the appearance of scatter in the N/O-O/H 
plane, according to the time-delay idea, namely that N is 
released some time after O. The time interval between the 
two bursts of a pair is tuned to give maximum scatter. It is 
found that this timescale is comparable to the timescale of 
star formation in an expanding supershell. The requirement 
for the 'shell burst' to be the second pair-burst is that the 
O-rich ejecta mix into the supershell. Because of the poorly 
understood physics of supershells and mixing processes, this 
assumption should be seen so far as a working hypothesis. 
All bursts are assumed to be instantaneous, hence repre- 
senting short but intense star formation events. The closed 
model is able to explain the observations of N/O-O/H in 
both scatter and level. If assuming the upper limit of hot 
bottom burning to be 5Mq, a=l.l is used to explain the 
observations, where a is the mixing length parameter. This 
is a rathe r low value compared to the results of recent works 



lations of stellar yields for metallicities lower than Z=0.008 
are definitely desired to quantify this.. 

A very important conclusion is that no primary N pro- 
duction in massive stars is needed to explain the observa- 
tions. Intermediate mass stars are in position to produce a 
sufficient amount of primary N. 

For the Y-O/H observations, a linear trend is visible, 
and a linear fit gives dY/dZ = 2.63 ± 2.21 and Y p = 
0.238 ± 0. 004. T hese values are consistent with those of Izo- 
tov et al. d!997a| ) within the uncertainties. The closed model 



how ever (van den Hoek & Groenewegen 1997; Marigo et al 
1996), Javouring a value close to 2. 

Using yield set 3 (the Padova set), it is necessary to 
extrapolate the primary N yields below Z=0.008 to obtain 
a higher N yield. Otherwise, the level of N/O becomes too 
low compared to the observations. The primary N produc- 
tion has to be increasing with decreasing metallicity. Calcu- 



is found to explain the Y-O/H observations perfectly, only 
with different slopes, depending on the yield set in use. In 
all cases the slopes are within the uncertainties. As seen b; 
inspection of the true He yields for the three sets in table 
one finds the explanation for the slope difference to be that 
the He yields of set 3 are 2-3 times higher than those of set 
2a or 2b. 

It is noteworthy that the Y-O/H relation was fitted us- 
ing exactly the same parameters as for the N/O-O/H fitting. 
The problem arises when fitting O/H-^i data. It is argued 
that the observed gas fractions are actually lower limits, be- 
cause dark matter is implicitly included in dynamical mass 
estimates and molecular hydrogen ignored. Hence, it is found 
useful to calculate upper limits using M/L estimates from 
starburst evolutionary models. For most objects, very ex- 
treme upper limits have to be used, assuming the galaxies to 
experience maximum luminosity of their first burst, except 
for three galaxies where the known star formation histories 
allow us to adopt more moderate and realistic upper limits. 
If star formation histories are found for a larger sample of 
dlrrs and even BCGs constraining chemical evolution mod- 
els with observed gas fraction intervals may eventually turn 
out to be extremely useful. The closed models do not repro- 
duce the gas fractions of the three well-known objects, even 
if a lower IMF-cutoff, equal to 0.01 M instead of 0.1 M 
is adopted. Hence, open models are considered, allowing gas 
to escape or to accrete on to the galaxy. Two kind of winds, 
enriched and ordinary, have been used. The results when in- 
corporating enriched winds are not in accordance with the 
observed level of N/O and the Y-O/H fitting is not satisfac- 
tory, when the gas fractions of the three well-known systems 
are fitted. Hence, it is concluded that our models employing 
enriched winds is in conflict with the observations. 

The next step is to include ordinary instead of enriched 
winds. Both N/O-O/H and Y-O/H are fitted, with results 
resembling those of the closed model. Unfortunately, only 
the gas fraction of one of the three selected galaxies is fit- 
ted, not differing much from the results of the closed model. 
To check the inclusion of ordinary winds, the outcome of 
the model is compared to the results of a simple analytical 
model, employing continuous star formation. A close resem- 
blance is found between the two models at low metallicities, 
but at higher metallicities, the numerical model seems to 
have problems in getting the yield down. The difference may 
be caused by the behavior of starbursts at a low absolute 
gas mass. This is confirmed when comparing the numeri- 
cal model, changed slightly to employ continuous star for- 
mation, to the analytical model, displaying almost identical 
outputs. Thus, it is concluded that it is important to spec- 
ify clearly whether instantaneous bursts or continuous star 
formation is used, when including ordinary winds. Instanta- 
neous bursts resemble the intense bursts of BCGs, whereas 
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the more moderate bursts of dlrrs are better explained using 
a continuous SFR. 

Finally, inflow and ordinary winds were included. The 
results are in accordance with all observations, except that it 
is difficult to obtain the right gas fractions and N/O scatter 
simultaneously. Only for quite extreme parameter choices 
as in fig. ^ and ^J, one may be succesful. ft may be im- 
portant to note that the results, when using the combined 
inflow/ordinary wind model, show the upturn in N/O to be 
more pronounced, than it was for the closed model. 

One question is unavoidable: which model is preferred? 
It is impossible to give an unambiguous answer. Dwarf galax- 
ies are different, both in mass and appearance. Some are ex- 
plained well by a closed model, others need a combination 
of ordinary winds and inflow. This is true for both dlrrs an 
BCGs. However, for dlrrs one should prefer to adopt a more 
continuous SFR before fitting the observations. 

In all cases, the model including enriched winds seems 
to be ruled out, since it is in direct conflic t with the obser- 
vations, as also found by e.g. Carigi et al. ( 1999 ). 
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